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ABSTRACT 



We consider a formation scenario for supramassive neutron stars (SMNSs) 
taking place through mass and angular momentum transfer from a close com- 
panion during a Low Mass X-ray Binary (LMXB) phase, with the ensuing sup- 
pression of the magnetic field. We show that this formation channel is likely to 
work for all equations of state except the stiff est ones, L, N*. After the end of the 
mass transfer phase, SMNSs will loose through magnetic dipole radiation most 
of their angular momentum, triggering the star's collapse to a black hole. We 
discuss the rate of occurrence of these collapses, and propose that these stars, 
because of the baryon-clear environment in which the implosion/explosion takes 
place, are the originators of gamma ray bursts. 

Subject headings: gamma rays: bursts - stars: neutron - black holes - relativity: 
general - instabilities 

1. Introduction 

One of the key requests on gamma ray burst (GRB) models is that they make contact 
with the fireball model (Rees and Meszaros 1992) which has proven so successful in 
predicting and interpreting the observed properties of GRBs' afterglows. In particular, this 
entails that a large explosion is to take place in a region with small baryon contamination: 
for E = 10 53 erg, the baryon contamination must be at most E j^c 2 fa 10 -4 M , for 7 = 300, 
the bulk Lorenz factor of the explosion. Vietri and Stella (1998) presented a model which 
could accomplish this, involving a supramassive neutron star (SMNS), i.e., a neutron star 
with a larger baryon number than any normal neutron star because it derives part of its 
support against self-gravity from the centrifugal force; these supramassive stars cannot 
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be slowed down to zero spin rate, because they are so massive that, as they lose angular 
momentum, they become unstable to black hole formation before reaching zero spin rate 
(Cook, Shapiro and Teukolsky 1994a,b, Salgado et al, 1994). The model consists of the 
implosion/explosion (I/E) of a supramassive neutron star which has lost through magnetic 
dipole radiation so much angular momentum that it must then collapse to a black hole; the 
rotational energy of the small amount of equatorial mass left behind because already in 
near centrifugal equilibrium provides the energy source that powers the burst. 

In paper I we proposed that SMNSs are formed in the SN explosion of a core with 
too much mass and angular momentum to end up in a normal neutron star. Though still 
nothing stands against this possibility, and we are not reneging it, we have now realized 
that a different channel exists: mass and angular momentum accretion from a companion 
in a Low Mass X-ray Binary (LMXB). The following discussion is relevant to the formation 
of MilliSecond Pulsars (MSPs), and will mimick arguments used in discussing the evolution 
of LMXBs into MSPs. We first discuss how the accretion of large amounts of mass and 
angular momentum may be realized in Nature and then apply this scenario to GRBs. 

2. Mass and angular momentum accretion 

The main obstacle to accretion onto a normal neutron star of large amounts of angular 
momentum from a close companion via an accretion disk is the neutron star's magnetic 
field: the neutron star can rotate only so fast as to make the corotation and Alfven radii 
coincide, lest a propeller phase sets in, which would actually entail angular momentum loss 
(Ghosh and Lamb 1978, Illarionov and Sunyaev 1975). The coincidence of these two radii 
leads to an equilibrium period, P eq = 1.3(5/10 12 G) 6 / 7 (M / M Edd )~ 3 / 7 s (Ghosh and Lamb 
1992), which clearly shows that B must decrease before significant spin-up can occur. 
Though no unique model has emerged yet, the current consensus is that the neutron star 
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magnetic field decays by at least 3 — 4 decades either as a direct result of mass accretion 
(Phinney and Kulkarni 1994) or of the ensuing spinup (Ruderman, Zhu and Chen 1998). 

The strongest constraints on field decay in NSs come from LMXBs and MSPs. In only 
one LMXB, SAX J1808. 4-3658 a coherent 2.5 ms signal has been detected in the persistent 
X-ray emission, providing direct evidence for the presence of a small magnetosphere; the 
inferred magnetic field is in the B ~ 10 8 — 10 9 G range (Psaltis and Chakrabarty 1998). 
All other LMXBs have undetectably small coherent pulsations in their persistent emission, 
if at all. Yet spin periods in the in the 2-4 ms range have been deduced for about ten 
LMXBs from the X-ray flux oscillations that are present during type I bursts emitted by 
these sources (cf. van der Klis 1998). Spinup through accretion can have occurred in these 
neutron stars only if their magnetic field is lower than ~ 10 9 G. Further evidence that the 
field might have decayed to dynamical insignificance derives from the modelling of the kHz 
Quasi Periodic Oscillations (QPO), a common phenomenon observed in LMXBs. In the 
sonic point model (Miller, Lamb and Psaltis 1998), the Alfven radius is located at a radius 
corresponding to a Keplerian frequency of ~ 350 Hz, corresponding to a magnetic field of 
w 8 x 10 8 G: this already bears witness to a thousand-fold reduction of the magnetic field 
below that of a typical newborn pulsar. A better model explains QPOs in terms of the 
fundamental frequencies of test particle motions in the general-relativistic potential well in 
the vicinity of the neutron star (Stella and Vietri 1999). The model is capable of explaining 
the observed relation between peak QPO frequency and its lower frequency counterpart over 
three orders of magnitude in peak QPO frequency and several distinct classes of sources, 
including candidate black holes and LMXBs (Stella, Vietri and Morsink 1999), implying 
that the magnetic field has been reduced already to < 2 x 10 8 G during the LMXB phase. 

A different argument involves a handful of MSPs with observed magnetic fields 
< 2 x 10 8 G , of which there are currently about a dozen, including the lowest fields ever 
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measured, 7 x 10 7 G in J2229 + 2643 and J2317 + 1439 (Camilo, Nice and Taylor 1996). 
For these small fields, the Alfven radius for disk accretion (Ghosh and Lamb 1978, 1992) 
is smaller than ~ 12 km, which, according to Cook, Shapiro and Teukolsky (1994b) is 
larger than the radius of the innermost stable orbit for the softest equations of state for a 
neutron star with M = 1.4 M Q (see Table 1). It should also be noticed that it has been 
argued cogently (Arons 1993) that the magnetic field in these objects is dominated by the 
dipole component, with negligible contribution from higher multipoles. Altogether, this 
means that we have already observed the result of accretion from a companion pushing the 
magnetic field to dynamical irrelevance (or, at least, very close to it) sometime during the 
mass exchange process. The argument about QPOs (Stella and Vietri 1999) implies that 
this may happen reasonably early in the LMXB history. 

Detailed models are required to establish the exact history of a neutron star's mass, 
angular momentum and magnetic field, but unfortunately these computations are currently 
fraught with uncertainties: where is the B field located, in the core or in the crust? And 
what is an appropriate model for the field suffocation? Population synthetic studies of 
this phenomenon (Possenti et al, 1999) have focused on two representative equations of 
state, and modeled the decay of the magnetic field in two limiting cases, imposing at 
the crust-core boundary either complete field expulsion by the superconducting core, or 
advection and freezing in a very highly conducting transition shell. The main result lies in 
the establishment of the existence of a tail in the rotation period distribution extending well 
beyond the shortest period observed so far (P = 1.558 ms), with only moderate dependence 
on the field suppression mechanism. For the softest equation of state the period distribution 
is still increasing at the shortest value before the onset of mass shedding, where Possenti 
et al. stopped their computations, while for the stiffest one the period distribution had a 
wide maximum around P = 2 — 4 ms, and a tail extending below this value. The fraction 
of objects with P < 1.558 ms is ~ 1% and rs 10% for the stiff and soft equation of state, 
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respectively, all ending up with very small magnetic fields, < 10 8 G. 

Though the accreted mass is larger when account is taken of the need to suppress 
the magnetic field than when the magnetic field is neglected, the difference is not very 
large (Burderi et al, 1999). So, in order to appraise whether the neutron stars thusly 
formed may be (or not) supramassive, we simply consider Table I, from Cook, Shapiro and 
Teukolsky (1994b). It shows that the total amount of mass that needs being accreted from 
a companion in order to reach the supramassive stage at the initial point of mass shedding 
depends strongly upon the equation of state. For equation of state C, the neutron star 
collapses to a black hole even before reaching mass shedding. The soft equations of state, 
EoSs, (A, D, E, KC) have become supramassive; the intermediate EoSs (C, M, UT, FPS) 
are within < 0.1 M Q of doing so, and will cross the threshold if accretion continues after 
the mass shedding point is reached (see below). The total amounts of mass required to 
become supramassive (~ 4/3 of the difference in mass at infinity, Phinney and Kulkarni 
1994, corresponding to a further ss 0.5 M Q ) are so modest that it seems likely that even the 
models based upon EoSs AU and UU will reach this stage, provided a donor of sufficient 
mass is found: this too will be discussed below. EoSs L and N* are hopeless: the total 
amount of mass to be accreted corresponds to ~ 2 M Q . If either of these EoSs were correct, 
there would be no way to form SMNSs via accretion from a companion in a LMXB. 

Both Cook, Shapiro and Teukolsky (1994b) and Possenti et al. (1999) halted their 
computations when the mass shedding rotation rate is reached, but there is nothing magic 
about this moment. Instead, Popham and Narayan (1991) and Paczynski (1991) argued 
that accretion continues unimpeded, in Newtonian stars plus disks configurations, with stars 
remaining close to the breakup angular speed, while mass and total angular momentum 
increase. Other reasonable possibilities may contribute to prolong mass accretion: the 
reduction of angular momentum through gravitational wave losses (propitiated by the 
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growth of a small stellar eccentricity) or the setup of a spiral shock wave reducing the 
angular momentum of incoming disk material, with the outward transport of angular 
momentum. In any case, several avenues are possible which would keep the neutron 
star marginally inside the mass shedding limit. For this reason, it seems nearly certain 
that the intermediate EoSs (C, M, UT, FPS) which are only ?s 0.1 M away from being 
supramassive, will reach this stage as mass accretion continues. 

Since all soft and intermediate equations of state only require ~ 0.5 M to become 
supramassive, their companion star may be any star with mass < 1 M , exactly as 
discussed in the normal recycling model for MSPs. EoSs AU and UT, instead, require 
in total ps 1 — 1.1 M to become supramassive. At first sight, this requirement might 
seem unsurmountable: thermal stability in the mass exchange process through Roche 
lobe overflow (Webbink et al. 1983) requires that the companion of the NS has a mass 
below 5/6 of the neutron star's, i.e., 1.17 M for an initial NS mass of 1.4 M . Since the 
smallest He-core that may be left behind is that of a star which took all Hubble time 
to evolve off the main sequence, 0.16 M , this leaves a maximum transferable mass of 
1.01 M , less than required for either AU or UT. However, this is incorrect: the famous 
requirement of 5/6ths only occurs because Webbink et al. (1983) considered Paczynski's 
(1967) approximation for the Roche lobe radius: using instead Eggleton's (1983) formula, 
this requirement disappears. Webbink et al. (1983, Eq. 15) show that the mass transfer 
rate is M x oc — I / (d In R^/d In Mi) where Rl is the Roche lobe radius, and they argue that 
thermal stability in the process requires X = din Rl/ din Mi > 0. Using Eggleton's formula 
Rl/o> — 0.49/(0.6 + ln(l + q 1 ^ 3 )/q 2 ^ 3 ) where a is the distance between the two stars and 
q = Mi/M 2 is the mass ratio, we find, under the hypothesis of conservative mass transfer, 

(1 + g)(2(l + q 1 ' 3 ) ln(l + q 1 ' 3 ) - q 1 ' 3 ) 

(l + g 1 /3)(l. 8g 2/3 + 31n ( 1 + g l/3)) <■ > 

for every q\ Also, dynamical stability exists provided the donor mass is < 2 M (Rappaport 



et al, 1995). So we may consider as a possible companion for neutron stars with EoSs 
AU and UT, sub/giants of mass < 2 M Q , where mass transfer is pushed forth by donor 
radius expansion, in complete analogy with the model of Webbink et al, 1983, except for 
donor mass. From Fig. 8b of Verbunt (1993), we see that a giant or subgiant of nearly 
solar metallicity, of, say, 1.7 M manages to transfer at sub-Eddington rates ~ 1.4 M to 
the neutron star, provided mass transfer begins when the giant core is ~ 0.2 M Q ; mass 
transfer will then leave behind a small (~ 0.3 M ), nearly inert He nucleus, with final 
period in the range of 0.3 d. We thus see that these systems provide attractive progenitors 
for supramassive neutron stars, even in the case in which the applicable EoS is either AU or 
UT, provided of course mass accretion is close to conservative, an implicit assumption we 
made throughout, and that mass can be accreted in sufficient quantities. 

Recent studies of binary pulsar masses (Thorsett and Chakrabarty 1999) seem to 
argue against significant mass accretion, but it should be noticed that, by investigating 
millisecond pulsars with periods exceeding ps 2 ms, the authors are investigating objects 
for which we know a priori that little mass need have been accreted, since their periods 
are long compared with SMNSs'. We may expect different results when pulsars are chosen 
otherwise: a recent redetermination of the mass of Cyg X-2 finds M = (1.8 ± 0.2) M & 
(Orosz and Kuulkers 1999), departing from the narrow range of Thorsett and Chakrabarty. 

The lowest magnetic field for the formation of supramassive neutron stars may be lower 
than the empyrical value (^ 2 x 10 8 G) mentioned above (Possenti et al. (1999) because 
when mass accretion from the companion begins to taper off, or alternatively if mass 
accretion is intermittent, the Alfven radius (which scales as M~ 2 / 7 ) may expand further 
than the corotation radius: the neutron star then goes through a new propeller phase which 
slows its rotation. The overall effect is not large so that we shall consider in the following a 
maximum magnetic field q x 10 8 G, with q < 1. 
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3. Further evolution 

Supramassive neutron stars are unstable to collapse to a black hole when angular 
momentum losses reduce the initial angular momentum to about half of the initial value; 
furthermore, these stars are peculiar in that evolution at constant baryon number, but 
decreasing total angular momentum, makes them spin up, rather than down; all of this 
is especially evident in Fig. 7-10-13-16 of Salgado et al, 1994. Magnetic dipole losses 
cause a net torque which spins down the neutron star in a time (Vietri and Stella 1998) 
tsd — 5 x 10 9 yr(10 8 G/B) 2 . This time-scale is not strongly dependent upon EoS, but 
depends strongly upon whether the model is only marginally supramassive, or close to the 
absolute maximum mass (rotating or not) for the given EoS, so that it may be considerably 
shorter under many circumstances. Thus, a time t s a after becoming supramassive, the 
neutron star will collapse to a black hole. This time is reasonably long when compared 
with typical mass accretion time-scales, which, as discussed above, are typically determined 
by sub/giant nuclear evolution timescales. Thus mass transfer will have long since ceased, 
and the immediate SMNS surroundings will be reasonably baryon free. The companion 
star, in the meantime, will have settled down as a low-luminosity, low-mass white dwarf, 
which is not expected to pollute the environment either. Furthermore, we can gauge the 
baryon-cleanliness of the SMNS surroundings at large if we assume that MSPs are born 
through the same chain of events, except less extreme, for then we know the Galactic 
distribution of MSPs. These are often located well outside the Galactic disk, within an 
ISM with typical densities well in defect of n — 1 cm -3 , which makes the total baryon mass 
within, say, 0.1 pc, less than 10 -5 M Q , more than enough to guarantee contact with the 
fireball model. We thus see that also this version of the formation scenario guarantees a 
baryon clean environment, exactly like the different scenario of Paper I. 

The situation is clean even in the case in which the collapse occurs while mass transfer 
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is still taking place. The total amount of baryons in the accretion disk is negligible: the disk 
crossing time is of order of « 1 month, which, with mass transfer rates 10~ 9 — 10~ 8 M Q yr^ 1 , 
corresponds to much less than the maximum contamination value. The total amount of 
outlying mass from a wind is also rather small: for M w fts 1O~ 9 M yr~ l and v w ~ 30 km s _1 , 
the total mass within, say, 0.1 pc is 3 x 1O~ 6 M , again negligible. The highly relativistic 
ejecta and 7 rays from the burst will hit the companion and form a shock way inside the 
star's photosphere, so that local dissipation of the ejecta kinetic energy will lead to the 
companion's inflating on the (long!) Kelvin-Helmholtz time-scale, and the non-thermal 
afterglow emission will not be contaminated by the re-radiated thermal component. 

The mechanism for the energy release is the same as discussed in paper I: once the 
neutron star is destabilized, the innermost regions will collapse promptly to a black hole, 
while the equatorial matter, which is close to centrifugal equilibrium, will just contract a 
little bit and begin orbiting the newly formed black hole. A necessary condition which needs 
to be met is that this equatorial material lies outside the innermost stable orbit. This can 
be checked from Table I of Cook, Shapiro and Teukolsky (1994b) who show that neutron 
stars which have reached the mass shedding regime have equatorial radii larger than the 
innermost stable orbit (see their column j), independent of EoS. In paper I, we estimated 
the amount of matter left behind as ~ 0.1 M Q ; this configuration is identical to the one 
hypothesized in most current models (Meszaros 1999), and the debris torus is massive 
enough to power any burst, especially in the presence of a moderate amount of beaming. 

We now discuss the rate at which spun-up SMNSs collapse to black holes. Since the 
timescales involved are a fair fraction of the age of the Universe, and since star formation 
evolves strongly in the recent past (Madau et al, 1996), we have to consider cosmological 
evolution of the population. However, from Fig.l of White and Ghosh (1998), it can be seen 
that the population of MSPs is roughly constant (within the accuracy of the present, order 
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of magnitude estimates) over the redshift range < z < 1 for most assumptions. There 
are currently an estimated 5 x 10 4 MSPs in the disk of the Galaxy (Lorimer 1995, Phinney 
and Kulkarni 1994); assuming that there are as many systems in the bulge, that a fraction 
/3 of these are SMNSs, and that the typical timescale for collapse to black hole is given 
by t s d, the expected rate of collapses in the Milky Way is r = 10 5 j3/t s d = (3/(5 x 10 4 yr), 
which is to be compared with the inferred rate of GRBs, 1 every 3 x 10 7 yr in an L± galaxy 
like the Milky Way. Scaling to (3 — 0.05, a value intermediate between the extremes of 
the simulations of Possenti et al, we find that the two rates agree for a beaming fraction 
(X1/47T ~ 0.2(/3/0.05); this is consistent with the idea that these explosions do not require 
extreme beaming fractions , since the explosion need not wade its way through a massive 
stellar envelope, but immediately breaks free into a baryon clean environment. 

This model makes an easily testable prediction, because the location of bursts inside 
their host galaxies is the same as that of LMXBs, which are distributed at distances from 
the Galactic plane z f« 1 kpc, most likely arising from kick velocities at the time of neutron 
star formation (van Paradijs and White 1995). A similar z distribution is observed for 
MSPs, z > 0.7 kpc, and moderate transverse speeds. Thus we would expect GRBs to 
cluster around galactic disks (contrary to the binary pulsar merger model, where at least 
some 50% of all GRBs should be uncorrelated with the original birth galaxies), but should 
not correlate with star forming regions (except for SMNSs which form directly during a SN 
event, as discussed in paper I), contrary to all scenarios involving massive stars. Also, the 
redshift distribution of GRBs within this model should be flatter than the star formation 
distribution (again contrary to hypernovae), because the redshift distribution of the MSP 
population is rather flat (White and Ghosh 1998). 

We acknowledge helpful discussions with G. Ghisellini, A. Possenti and L. Burderi. 
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Accretion from Keplerian disk onto a 1.4 M neutron star 



EoS a 


M f /M b Q 


Mi/M C Q 


AM/M| 


pe 


R-i 


K 


M sm jM% 


A 


1.77 


1.57 


0.428 


0.604 


9.59 


2.82 


1.66 


C 


1.74 


1.54 


0.389 


0.894 


12.1 


0.27 


1.86 


D 


1.76 


1.56 


0.405 


0.730 


10.7 


1.70 


1.65 


E 


1.76 


1.57 


0.414 


0.656 


10.0 


2.38 


1.75 


F 


1.52 


1.59 


0.172 


0.715 


9.21 


3.20 


1.46 


L 


1.80 


1.52 


0.443 


1.250 


15.0 


0.00 


2.70 


M 


1.74 


1.50 


0.367 


1.490 


16.7 


0.00 


1.80 


N* 


1.84 


1.53 


0.484 


1.080 


13.6 


0.00 


2.64 


KC 


1.74 


1.55 


0.385 


0.888 


12.1 


0.31 


1.49 


AU 


1.79 


1.58 


0.446 


0.701 


10.4 


2.01 


2.13 


uu 


1.78 


1.56 


0.436 


0.784 


11.2 


1.26 


2.20 


UT 


1.78 


1.57 


0.429 


0.754 


10.9 


1.52 


1.84 


FPS 


1.76 


1.56 


0.416 


0.747 


10.9 


1.56 


1.80 



Table 1: From Cook, Shapiro and Teukolsky, 1994b. a Equation of state; b final total mass- 
energy: c initial rest mass; d accreted rest mass; 6 rotation period in ms; * initial circumfer- 
ential radius in km; 9 initial circumferential height of corotating marginally stable orbit in 
km; h maximum static total mass-energy for EoS. 



